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Abstract 

Graphical models of probabilistic dependencies have been 
extensively investigated in the context of classical uncer- 
tainty. However, in some domains (most notably, in com- 
putational physics and quantum computing) the nature of the 
relevant uncertainty is non-classical, and the laws of classi- 
cal probability theory are superseded by those of quantum 
mechanics. In this paper we introduce Markovian Entan- 
glement Networks (MEN), a novel class of graphical repre- 
sentations of quantum-mechanical dependencies in the con- 
text of such non-classical systems. MEN are the quantum- 
mechanical analogue of Markovian Networks, a family of 
undirected graphical representations which, in the classical 
domain, exploit a notion of conditional independence among 
subsystems. 

After defining a notion of conditional independence appropri- 
ate to our domain (conditional separability), we prove that the 
conditional separabilities induced by a quantum-mechanical 
wave function are effectively reflected in the graphical struc- 
ture of MEN. Specifically, we show that for any wave func- 
tion there exists a MEN which is a perfect map of its condi- 
tional separabilities. Next, we show how the graphical struc- 
ture of MEN can be used to effectively classify the pure states 
of three-qubit systems. We also demonstrate that, in large 
systems, exploiting conditional independencies may dramat- 
ically reduce the computational burden of various inference 
tasks. In principle, the graph-theoretic representation of con- 
ditional independencies afforded by MEN may not only fa- 
cilitate the classical simulation of quantum systems, but also 
provide a guide to the efficient design and complexity analy- 
sis of quantum algorithms and circuits. 



Introduction 

Probabilistic inference in large systems is a topic of great in- 
terest in computer science. In contexts where the underlying 
uncertainty is of classical nature, the graphical representa- 
tion of probabilistic dependencies and independencies has 
proved to be a valuable tool in order to reduce the computa- 
tional burden of various inference tasks (Pearl 1988). Two 
main classes of graphical representations are of special in- 
terest in the classical domain. Bayesian Networks (BN) are 
based on directed, acyclic graphs which capture causal de- 
pendencies among variables of interest. By contrast, Marko- 



vian Network^] (MN) are non-causal representations based 
on undirected graphs, which encode mutual symmetric de- 
pendencies among the variables of interest. Whereas the BN 
formalism is especially convenient at the modeling stage, as 
causal dependencies among variables of interest are more 
easily identified, the structure of MN is particularly advan- 
tageous from a computational viewpoint. In particular, sev- 
eral algorithms for probabilistic inference in BN require as 
a first step the transformation of the original BN into a cor- 
responding MN (e.g., the junction tree algorithm). 

Recently, there has been some interest in the graphical 
representation of non-classical uncertain domains, most no- 
tably in the context of computational physics and quantum 
computing. (Tucci 1997) introduced Quantum Bayesian 
Networks (QBN), which provide a structured representation 
of quantum-mechanical wave functions along the lines of 
classical BN. 

Missing a quantum-mechanical generalization of MN, the 
algorithmic toolbox of QBN is significantly more limited 
than that of BN. Also motivated by such limitations, we in- 
troduce a novel class of undirected networks, Markovian En- 
tanglement Networks (MEN), which extends the formalism 
and computational gains of MN to quantum-mechanical sys- 
tems. 

We also introduce a novel notion of conditional indepen- 
dence appropriate to the context of entangled systems (con- 
ditional separability), which generalizes the classical notion 
of conditional probabilistic independence, and provides a 
rigorous semantics for our graphical framework. 

While in applications which only involve classical uncer- 
tainty the computational performance of MEN is equivalent 
to that of MN, many of the computational benefits of the 
modular structure of MN carry over to MEN even in appli- 
cations characterized by non-classical uncertainty. In par- 
ticular, the graphical structure of MEN can in principle be 
exploited in the context of classical simulations of quantum- 
mechanical systems (e.g., protein folding), and as a guide 



'The word Markovian in graphical representations does not re- 
late directly to Markov processes (a family of stochastic processes 
in which the probability distribution of future states for a given 
history only depends upon the current state, i.e., future states are 
conditionally independent of past states given the present state) but 
rather to the Markovian property of conditional independence of a 
random variable from all other variables given its neighbors. 



to the efficient design and complexity analysis of quantum 
algorithms and circuits. 

The paper is organized as follows. In the second section 
we introduce conditional separability, and relate it with both 
separability and the classical notion of conditional proba- 
bilistic independence. In the third section we formally intro- 
duce MEN, and discuss some of their structural properties. 
In particular, we show that the conditional separabilities in- 
duced by any wave function are precisely characterized by 
the graphical structure of an appropriately selected MEN. 
In the fourth section we address the issue of probabilistic 
inference, showing how marginal and conditional probabili- 
ties can be recovered from the structural elements of MEN, 
and discussing the computational performance of MEN in 
the context of some typical inference tasks. The fifth sec- 
tion provides a classification of 3-qubit states based on their 
network topologies, while the last section concludes. 

Setting, a- Independence, Separability, 
Conditional Separability 

Let us first introduce the notion of a quantum bit (qubit), 
the fundamental unit of quantum information, which can be 
thought of as a quantum mechanical analogue of a classical 
data bit. A (pure) state of a quantum bit is represented by 
a vector in a two-dimensional complex Hilbert space. The 
two computational basis states are conventionally written as 
|0) and |1) , and form an orthonormal basifl for the Hilbert 
space. A qubit state is a superposition of those basis states: 

|^) = a(|0))|0)+a(|l))|l), 

where a(|0)) and a(|l)) are complex numbers denoting 
probability amplitudes in the directions given by |0) and 
respectively. These amplitudes are normalized so that their 
Euclidean norm is unitary, i.e., |a(|0))|'^ + |a(|l))|^ = l.The 
reason is that the probabilities that the qubit will be observed 
in states |0) and |1) , respectively, are represented by the 
square moduH |a(|0))r and |a(|l))r, so with normalized 
amplitudes the total probability that the qubit is observed to 
be in either state |0) or |1) is 1. For notational convenience, 
in the remainder of this paper we shall give up the bra-ket 
notation |-) and always write x in place of\x). 

A pure state (wave function) of an n-qubit system takes the 
following form: 



V'. 



N 



a(x)x, 



where: 



N — { 1 , . . . , 7i} is a set of ?i qubits. 



^An orthonormal basis of a vector space V equipped with inner 
product (■, ■) is a subset {vi, Vk} of V which spans the whole 
space, and is such that (Vi, j) [{vi , Vj ) 
necker Delta. 



• X — xi ® ... (E) a;„ is a vector in an orthonormal basis 
X of the n-qubit system, while each Xi G {Oi, 1^} is an 
orthonormal basis of the i — th qubit, i = I, ...,n, and0 

• a{x) is a complex number denoting the probability am- 
plitude of the n-qubit system in the direction given by x, 
normalized so that J^xgx = 1- 

In this paper we only deal with pure states, while deferring 
the case of mixed statefl to future work. 

The state of a composite system is said to be separable 
(or not entangled) if it can be written as tensor product of 
states of the component systems (Nielsen and Chuang 2000). 
Formally, 

i'N = -ipM.N-M = iJ-'M 'E ll^N~M- 

Intuitively, the subsystems do not interact with each other, 
so they can be considered separately. This implies proba- 
bilistic independence between subsets of random variables 
across such non-interacting subsystems. 

Let {Xi}^^j^ be afinite, ordered set of Boolean random vari- 
ables representing the possible realizations of the n qubits, 
and let x^'^ = {x^, where 2;^° is a given realiza- 

tion of the i — th random variable, be an arbitrary reference 
pointQ. A joint realization a: AT = {xi,...,Xn) corresponds 
to a basis state of the composite n— qubit system. For any 
M C N we denote by a7 the set N - M, and by Xm the 
set {^ijjgM ■ 

Finally, we define a notion of a— independence which, as 
we shall see, is equivalent to separability. 

Two sets of qubits, M and M, are said to be a-independent 
if the following condition is satisfied for all xm and xjj: 

a(xM , xjj)a{x%,x^) = a{x%j,xji)a{xM,x^). 



Theorem 1 Two subsystems M and M are a-independent 
if and only if they are separable. Formally, 

V(a;Ar) a(xM,xjj)a(x%,x^) = a{x% , xjj)a(x m , x^) 

Proof. We prove Theorem 1 through a two-step proce- 
dure. We first show that separability of the two subsystems 
implies their a-independence. Next we show the converse, 
namely that a-independence implies separability. The two 
steps are sufficient to conclude the proof. 

Step 1. separability ^ a-independence 

If the two subsystems, M and Af, are separable then 



5ij], where 5ij is the Kro- 



The notation v®u stands for the tensor product of two vectors 
V and u. 

"^A mixed state is a classical probability distribution over pure 
states. Mixed states are conventionally represented by a density 
matrix p (i.e., a self-adjoint, positive-semidefinite matrix of trace 
one) with rr(p^) < L 

'We always use capital letters to denote random variables, and 
small letters to denote their realizations. 



i>N = i> M ® i^M- Let us write each of these wave functions 
in vector form: 



\ «(4 ') / 



where is the cardinality of set F, and 2l^l is the num- 
ber of basis states of the |F|— qubit system. 
Next, we rewrite the separabihty condition in the following 
form 



O-K^M ) 



V «(4 ) } 



or, equivalently. 



I 1 2l"U 



M 



M 



It is straightforward to verify that the following condition 
must then hold: 

for all x%[ , x\^ , x^ and x-!^. In particular, the condition 

must hold true when x^^ = x^ and x^ = Xj^. It foUows 
that separabihty implies a-independence. 

Step 2. a-independence^ separability 

By a— independence of systems M and M, we know 
that the following condition holds true for all xm and xjj: 

This can then be rewritten as 



a{x''}^,xL-) — 



M' 



a{x%) 



where a{x^) stands for the amplitude of the arbitrary refer- 
ence point x%, which without loss of generality we assume 
to be non-zero. 

Thus, the state of the n— qubit system takes the form 

/ a{x\j,xl^) \ 



JV 



a{x\^,xjj ) 



/ a{x\j , Xjj)a{x1[ , x^) \ 



1 



a{x\f , x^)a{x1f , x^ ') 



V aixr--ya{^,xll^'^) ) 
which factorizes in the following way 

"■\Pn) \ , 2lMI \ 



Let us denote the vectors in the above tensor product as a 
and /3, and rewrite the last equation as follows: 



1 



N 



a{x%) 



a® (3. 



Let Cq , Cf3 be complex numbers such that CaCf) = ^^^n-y, and 
\ca\ \\a\\ = \c(}\ ||/?||. It follows that 

where a' = c^a, /3' = Taking the norm on both sides, 
and using the fact that ipN has unit norm, we find that 



ll"'ll 



1. 



This imphes ||q:'|| = = 1 (i.e., a', and /?' are pure 
states of the component sy stems). ■ 

The above characterization of separability via 
a— independence makes it possible to introduce in a 
natural way a notion of conditional separability, defined as 
follows. Let {A, B, C'} be a partition of N into non-empty 
subsets. We say that A and B are conditionally separable 
given C if 

\/{xN)[a{xA ,XB,xc)a{x%x%,xc) 

= a{x% XB,xc)a{xA, x%,xc)]- 

Whenever two subsystems are not conditionally separable, 
we say that they are conditionally entangled. 

Theorem 2 Conditional separability implies conditional 
probabilistic independence. 

Proof. Observe that conditional probabilistic indepen- 
dence of Xa and Xb given Xc can be expressed as 



V(xjv) [p{xa\xb,xc) = p{xa\x%,xc) 



where 



p{xa\xb,xc) 



p{xa,xb,xc) 
T,^.p{xa,xb,xc)' 



In terms of probabihty amphtudes, this condition takes the 
following form: 



By multiplying the above fraction's numerator and denom- 
inator by | a (x^^, x^, xc) | and by the fact that I ci II C2 1 = 
IC1C2I, where ci, C2 G C, we obtain 



p{xa\xb,xc) 



\a{xA,XB, xc)a{x°J^,x%, xc)\ 



Y.XA \a'{xA,XB,xc)a{x\,x%,xc)\ 



2 ■ 



Next, from conditional separability of A and B given C, one 
obtains that 



p{xa\xb,xc) 



\a{x\, XB,xc)a{xA,x%,xc)\'^ 
ExA \a{x%XB,xc)a{xA,x%,xc)\'^ 

^ \a{x\,XB-,xc)\^\a{xA,x°B,xc)\'^ 
|a(a;^,a;B,a;c)PExA \a{xA,x%,xc)\'^ 
^ \a{xA,x%,xc)\'^ 
T,XA \a{xA,x%,xc)\'^' 

Observing that the latter term is equal to p{xA\xg , xc) con- 
cludes the proof. ■ 

Finally, we introduce a potential function that we shall uti- 
lize in defining Markovian Entanglement Networks. Follow- 
ing the approach of (La Mura and Shoham 1999), we define 
the amplitude potential function {q— function) as follows: 



q{xm\xm) 



a{xM, x^) 
a{x\j,xjj) 



The g-function, whenever it is well defined, can be inter- 
preted in terms of ceteris paribus comparisons: it describes 
how the probability amphtude changes when the reahzations 
of the qubits in _M are shifted away from the reference point, 
while those in M are held unchanged at xjj. 

Note that, whenever q is well defined, the condition 

\/{xN)[q{xA\xB,xc) = q{xA\x°BTXc)\ exactly corre- 
sponds to conditional separability of A and B given C, 
in which case we shall simply write q{xA\xc) in place of 

q{XA\XB,Xc). 

Markovian Entanglement Networks: a formal 
definition 

We define a Markovian Entanglement Network as an undi- 
rected graph G{V, E) with nodes i E V representing 
quantum-mechanical subsystems, and edges {i, j} G E rep- 
resenting conditional entanglement between subsystems i 
and j. For simplicity of exposition, we proceed under the 
assumption that each node represents a single quantum bit. 
Finally, each node is associated with a non-zero function 
q{xi |a:j/(j) ) (as defined in the previous section), where U {i) 
denotes the set of nodes directly connected to i via entangle- 
ment edges (or neighbors). Figure 1 depicts a simple MEN. 

If the g— functions are specified directly, then any ar- 
bitrary assignment of non-zero complex-valued functions 
q{xi\xu_ (^i-) , x'^^^-^) for all i (where U+{i) denotes the set 
of all qubits in U{i) whose index is greater than i, and 




q(x2|x„ x/, X4«) 



qr(xjxj,. X|) 



Figure 1 : A simple MEN. Nodes represent qubits, and edges 
the fact that two qubits are conditionally entangled. 



U-{i) = U{i) — [/+(«)) uniquely identifies a correspond- 
ing wave function, up to the phase of the reference point. In 
fact, once all the g— functions are specified, the relative am- 
phtudes a{xN)/cL{x^j^) are determined as follows: 

a{xN) a(xi ,x^) a(xjv) 

~ a(x^) a(xi,x'^) 



a(xiv) 



= liXl (,) )g(x2 |X1 , 4^ H)) aixtZl^,) 

= ... = ll^q{x^\xu_(^),x'^J^^^^), 

where xrjXr' denote the realizations of the remaining 
qubits in the appropriate context. 

Since the state of the composite n-qubit system is pure, we 
obtain the following condition: 



orm 

\a{x%) 



Y{i(l{xi\xu_(i),x^u^^i))a{x%) =1. 

Performing some elementary algebra, we conclude that 

1 



Exjv Y{i(l{x^\xu.(i),xl^^^^) 



From the last condition follows that only the the modulus, 

but not the phase, of the reference amplitude a{x^j^) is iden- 
tified, which means that the corresponding wave function 
is only determined up to the phase of the reference point. 
Conversely, any wave function in a given computational ba- 
sis can be represented by an appropriately selected MEN, as 
long as all the amplitudes are non-zero. 
We remark that, if q{xi\xi^^s^i^) only depends on xui^i) , then 
fixing xif(i-j completely specifies the behavior of the wave 
function along the i — th coordinate (up to the phase of the 
reference point), and that such behavior does not depend on 
the specific values taken by the remaining variables. 
It turns out that node separation with respect to MEN charac- 
terizes all the impUed conditional separabilities. More pre- 
cisely, for any wave function there exists an undirected graph 
G such that, for any partition of N into three non-empty sets 



of quantum bits A, B, and C, A is conditionally separa- 
ble from B given C if, and only if, C separates A from B 
in graph G, i.e., every path from a node in A to a node in B 
passes through some node in C. In the terminology of (Pearl 
and Paz 1987), such a graph is said to be a perfect map of 
the independence structure. 

Theorem 3 The set of conditional separabilities generated 
by any wave function has a perfect map. 

Proof. We appeal to a necessary and sufficient condition 
in (Pearl and Paz 1987). Specifically, we check whether 
conditional separability satisfies the following five proper- 
ties: symmetry, decomposition, intersection, strong union 
and transitivity. 

Let A, B, C, D, R, R', R" be subsets of qubits, where 
R, R', R" denote the subset of remaining qubits in the ap- 
propriate context. 

For the purpose of this proof, let us say that A is condition- 
ally independent of B given C, and write /(A, B\C) if and 
only if 

iyxN)[a{xA:XB,xc,xji)a{x\, x%,xc, xr) 

= a{x\,XB,xc,XR)a{xA,x%,xc, xr)]. 

Observe that this formulation reduces to conditional separa- 
bility if A, B, C are a partition of N. 
Then the following properties hold. 

Symmetry: I{A,B\C) ^ I{B,A\C). 

This follows trivially from the definition of conditional 
independence. Specifically, I{A, B\C) ^ I{B, A\C). 

Decomposition: IiA,B U D\C) =^ I{A,B\C) A 
I{A,D\C). 

This is equivalent to 

[a{xA, XB,Xc,XD,XR)a{x\,X%,Xc, X^jj,Xr) = 

a{x°^,XB,xc, XD,XR)a{xA,x%,xc,x°jj,XR)] ^ 

[a{xA ,XB,Xc, Xr- )a{x\ ,Xb,Xc,Xr-) = 

a(x5i ,XB,xc, xr> )a{xA,x%,xc, xr/ )] A 

[a{xA, Xc, XD, XR'')a{xA,Xc, X%,Xr-) = 

a{x\ ,xc,xd, XR/')a{xA ,xc,x'jj, xr/- )] . 
This follows trivially, since R' = {D, K) and R" = {B, R). 

Intersection: 

liA, B\C U Z?) A I{A, D\B U C) ^ I(A, B U D\C). 
Equivalently, 

[a{xA,XB,xc,XD,XR)a{x"j^,x%,xc, xd,xr) = 

a{x\ ,xb,xc,xd, XR)a{xA ,x%,xc ,xd,xr)\ A 

[a{xA:XB,Xc,XD,XR)a{x\,XB:Xc, X°jj,Xr) = 

a{x°^,XB,xc,XD,XR)aixA, xb,xc, x°jj,xr)] ^ 
[a{xA,XB,xc,XD,XR)a{x°j^,x%,xc, x°jj,xr) = 



a{x^j^,XB,xc,XD,XR)a{xA,x%,xc,x%,XR)]. 
This follows quite easily by algebraic manipulation. 

Strong union: I{A, B\C) => I{B, A\C U D). 
Equivalently, 

[a{xA,XB,xc,XR)a[x%x%,xc,XR) = 

a{x'A,XB,xc,XR)a[xA,x%,xc,XR)] ^ 
[a{xA, xb,xc,xd, xw)a{x"j^ ,x%,xc,xd,xr>) = 

a{xA ,Xb,Xc,Xd, Xri )a{xA:X%,Xc,XD,XR>)\. 

This follows by symmetry, and the fact that R= {D,R'). 

Transitivity: I{A, B\C) => I{A, V\C) V /(B, V\C) 
where V is one qubit. This is equivalent to 

[a{xA,XB,xc,XR)a{x°^,x%,xc,XR) = 

a{x°j^,XB,xc,XR)a{xA,x°B,xc,XR)] ^ 

[a{xA ,xc,xv, XR')a{x°j^ ,xc,Xv, xr> ) = 

a{xA ,xc,xv, XR' )a{xA ,xc,Xy, xr> )] V 

[aixB,xc,xv, XR>,)a{x% ,xc,Xy, xru ) = 

a{x% ,xc,xv, XR" )a{xB,xc,Xy, xr" )] . 



which follows by observing that V is either an element of 
A, or B, or R. Thus, we appeal to Pearl and Paz's result to 
conclude the proof. ■ 

The relevance of this result lies in the fact that one can repre- 
sent and reason about conditional separabilities graphically 
(i.e., one can exploit the graphical structure of the perfect 
map to reason about conditional separabilities). 

Exact Probabilistic Inference in MEN 

Markovian Entanglement Networks contain information 
about quantum probabilistic dependencies, but the informa- 
tion is not explicit. It is encoded in the potential functions, 
together with the topological structure of the network. The 
basic operation of computing marginal probabilities, and 
hence conditional probabilities, can be reduced to manip- 
ulations of the wave potentials. In fact, one readily obtains 
p{xm)/p{x%) (where p{xm) is the marginal probability of 
obtaining xm) by summing up the squared moduli of the 
g-functions: 



Ut9{Xt\x 



The computational performance of MEN in evaluating 
probabilistic queries is similar to that of MN. In particular, 
as in the case of MN, exact inference in singly-connected 
network^ can be performed in time linear in the number of 
nodes. We provide time-complexity analyses for some basic 
inference tasks in the simplest singly-connected networks, 
chains (Figure 2). 




Figure 2: A MEN chain. 



= argmaxx^ . ..argmax^^ \q{xi\x^) ■ • • g(a;„|a;„_i)| 

= argmaxxi \q{xi\x2)\^ ■ ■ ■ argmaxx^ |g(x„|x„_i)|^ 

Thanks to the exploitation of conditional separabihties, the 
problem enjoys a decomposition which, once again, allows 
to find a solution in time linear in the number of qubits. 
Specifically, a global maximum can be found by solving a 
sequence of n simpler maximization problems. 



Specifically, we consider three types of queries: (i) com- 
puting marginal probabihties, {ii) computing conditional 
probabiUties given evidence, and {Hi) finding a maximum 
Ukelihood instantiation. 

(i) Computing marginal probabiUty ratios takes the 
following form: 



p{xm) 



E 



]j5(xi|a;i_i,a;°+i) 



The right-hand side of the above formula can be rearranged 
as follows: 

J2xN-M ni=l q{^i\Xi-lj ^i+l)q{Xm+l\Xm, ' ' ' 

■ ■ ■ q{Xn\Xn-l)\'^ 

We assume, for simplicity, that M C A'' is the set of qubits 
with indices from 1 to m, whereas in — M are all the 
qubits with higher indices. Hence, the following holds: 

Ex„+i ■•■Ex„ \UiLll(^i\^i--i-'^i+l) ■ ■ ■ <ii^n\x„-l)\^ = 

\^i-lT^i+l)\'^'^Xm+i |Q'(^™+ll^™'^m+2)| ■■■ 
■■■T,xn \q{Xn\Xn-l)f ■ 

The above decomposition makes it possible to com- 
pute the marginal probability of interest performing 
6(n — m) + 2m — 1 operations. Thus, the time complexity 
of this task is hnear in the number of qubits n. 

{a) Computing conditional probabilities given some 
evidence xm takes the form 

P[Xn-M\Xm) - ^(^77) - p^xm)Mx"^) • 

Since the computation of marginal probability ratios has lin- 
ear time complexity, this is also the case for the computation 
of conditional probabilities. 

(iii) Finding a maximum likelihood instantiation takes the 
form 



argmaxxfj \a {xn)\ argmaxx^ 



]Jg(xi|a;i_i,x°+i) 



A network is said to be singly connected if, in the associated 
graph, there exists at most one path between any two nodes. 



MEN for 3-Qubit Systems 

From Theorem 3, we know that all conditional separabilities 
induced by any given wave function can be encoded in the 
graphical layer of an appropriately chosen MEN. In the case 
of a two-qubit system, the dependencies generated are triv- 
ial: either the two qubits are entangled, or they are separa- 
ble. Hence, for two-qubit systems there are only two classes 
of MEN topologies.- one with cormected, and one with un- 
connected nodes. Since in two-qubit systems conditional 
separability plays no role, in this section we investigate con- 
ditional separabilities in the context of the simplest interest- 
ing class: three-qubit systems. First, following (Duer, Vi- 
dal, and Cirac 2006), we classify states of such systems into 
equivalence classes under orthonormal basis changes. Next, 
we analyze the topological properties of the corresponding 
MEN, and what structural changes they undergo once we in- 
stantiate (measure) a component qubit. 
(Duer, Vidal, and Cirac 2006) classifies three-qubit system 
states into fully entangled (i.e., states of composite systems, 
in which no component system is separable) and not fully 
entangled. Each of these classes is comprised of two sub- 
classes. The first one contains GHZ-like and W^-Uke states, 
while the second contains completely separable (i.e., there 
is no entanglement present) and bi-separable states (i.e., two 
qubits are entangled and separable from the remaining one). 
The classification is robust with respect to local operations: 
in particular, by means of basis changes no state in one class 
can be obtained from a state in a different class. It is a natu- 
ral question whether these four classes of three-qubit states 
can be characterized purely in terms of their MEN topology. 
All possible MEN topologies corresponding to the four 
classes of three-qubit states are summarized in Figure 3. 
One can easily see that the MEN topology of all classes ex- 
cept for GHZ-hke is invariant under basis changes. How- 
ever, in the case of GHZ-like states, there exist bases in 
which the MEN graph assumes each of the configurations in 
Figure 3(d) (i.e. the complete graph, or any permutation of 
a three-node chain). In passing, we observe that conditional 
separabihties, unlike unconditional ones, are not preserved 
under basis changes. We conclude that the four classes of 
three-qubit system stales inlroduced above are completely 
characterized by the topology of their associated MEN. 

Next, we consider the issue of measurement and its im- 
pact on MEN topology. Measuring (instantiating) a qubit 
collapses the wave function along one of the qubit' s ba- 
sis states, destroying any entanglement with the remaining 
qubits. This has an impact on MEN topology: specifically, 
after measuring a qubit, the MEN graph in the new state 



a) Separable state b) Bi-separable state 
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Figure 3: MEN for all classes of three-qubit states. 

is obtained by erasing all the edges which involve the mea- 
sured qubit. While measurement may destroy edges, it never 
introduces new ones, as conditional separabilities are always 
preserved under qubit instantiation. 

Conclusions 

We introduced Markovian Entanglement Networks, a novel 
class of graphical representations for quantum-mechanical 
systems, and argued that the modular representation of 
quantum-mechanical states gives rise to computational ad- 
vantages analogous to those afforded by classical Markovian 
Networks. We also showed that, based on MEN topology 
alone, one can effectively classify all possible three-qubit 
system states. The most related approaches in the literature 
are (Tucci 1997) and (Buzek and Plesch 2003). (Tucci 1997) 
introduced Quantum Bayesian Networks (QBN), a graph- 
ical representation based on directed acyclic graphs. In 
QBN the emphasis is on modeling, rather than on computa- 
tion, and QBN cannot be directly used to facilitate computa- 
tional tasks. (Buzek and Plesch 2003) introduced Entangled 
Graphs (EG), which capture the distribution of bipartite en- 
tanglement in multi-qubit systems. The graphical decompo- 
sition in EG is based on concurrence (a statistical measure of 
bipartite entanglement), and cannot directly represent condi- 
tional separabilities. Furthermore, missing a modular repre- 
sentation of amplitudes, EG cannot encode all the informa- 
tion contained in a quantum-mechanical state. 
Two main areas of application for MEN come to mind: 
the classical simulation of quantum-mechanical systems, 
and quantum computing. In particular, we hope that the 



identification and exploitation of conditional separabilities 
in quantum-mechanical systems may lead to the design of 
more efficient classical simulations in a variety of compu- 
tationally demanding applications, from physical chemistry 
to quantum optics. Furthermore, in the context of quantum 
computation, we hope that the graphical structure of MEN 
may provide a guide for the efficient design and complexity 
analysis of quantum algorithms and circuits. 
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